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Experimental results obtained with the HERA collider and recent progress in their theoretical interpretation 
are reviewed. After a short introduction to HERA physics, deep inelastic scattering and photoproduction are 
discussed as (virtual) photon-proton scattering. It is shown that the measurement and theoretical understanding 
of both photoproduction as well as low- a; deep inelastic scattering are essential for obtaining reliable high energy 
extrapolations within hadron-hadron interaction models. Limitations of the predictive power of hadron interaction 
models due to the interplay of perturbative QCD and unitarity effects are discussed. 



1. Introduction 

It is very difficult to compare data measured 
at the HERA collider directly to data obtained 
in cosmic ray experiments. Almost all cos- 
mic ray experiments deal with final state par- 
ticles (hadrons, electrons, photons, muons, or 
neutrinos) which are produced in interactions 
of primary cosmic ray (CR) particles with air. 
Furthermore, the experiments are very sensitive 
to the projectile fragmentation region. In the 
contrary, modern collider experiments focus on 
particle production in the central pseudorapid- 
ity range in the projectile-target center-of-mass 
system (CMS). They either deal with leptons, 
(anti)protons or heavy nuclei as primaries. Never- 
theless, it is possible to relate particle production 
processes in CR interactions to those studied in 
collider experiments. 

Since the number of final state particles is of- 
ten very large, efficient approximations as well as 
simulation techniques have to be applied to de- 
rive conclusions from measured data or to study 
theoretical predictions. One of the most powerful 
technique is the realization of theoretical mod- 
els in Monte Carlo event generators. Such event 
generators allow us to study hadron production 
at colliders as well as in CR interactions. 
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Unfortunately, our current understanding of 
hadronic interaction processes is rather limited. 
Therefore, these event generators can only be 
built by combining theoretical predictions with 
phcnomcnological models and parametrizations 
and have to be tuned by comparing their predic- 
tions to fixed target and collider measurements. 
Considering the underlying theory and models 
entering such Monte Carlo programs, almost all 
data measured at HERA are relevant to our un- 
derstanding of very high energy cosmic ray inter- 
actions. Examples include data on 

(i) limits on physics beyond the Standard Model, 

(ii) parton densities, 

(iii) reliability and range of applicability of 
perturbative QCD and various approximation 
schema, 

(iv) transition between soft and hard physics 
(perturbative and non-perturbative regimes), 

(v) forward hadron production, and 

(vi) heavy flavour production. 

It is beyond the scope of this work to discuss all 
these data. Instead, we shall focus on the impli- 
cations of the low- a; structure function measure- 
ments and recent results on forward jet |||| 
and leading baryon production . Where pos- 
sible, comparisons to predictions of Monte Carlo 
models currently used in the analysis of CR data 
are included. 

After a short introduction to the HERA col- 
lider, deep inelastic scattering (DIS) and parton 
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densities at low x are discussed in Sec. In 
Sec. U DIS is interpreted as (virtual) photon- 
proton scattering and data on leading jet and 
baryon production are presented. Finally, in 
Sec. U a summary is given. 

2. The HERA collider 

The electron-proton storage ring HERA at 
DESY, Hamburg (Germany) went into operation 
in 1992. At HERA, 27.6 GeV positrons or elec- 
trons are collided on 820 GeV protons. The lep- 
tonQ-proton cm. energy is \fs ep ~ 300 GeV which 
corresponds to a lepton energy of £7i a b w 47 x 10 12 
eV in the proton rest frame. 

The typical scattering process studied at 
HERA is sketched in Fig. [jj The beam lepton 




There are only two independent variables, we 
shall use x and Q 2 . In the proton infinite- 
momentum frame, x denotes the momentum frac- 
tion of the proton carried by the struck quark. 

In Fig. U the kinematic range for DIS accessible 
by various fixed-target experiments and HERA 
is shown. The HERA measurements extend the 
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Figure 2. Kinematic region of DIS covered by 
experiments up to now Q. 



Figure 1. Kinematics of positron-proton scatter- 
ing in Breit frame. 



scatters off the proton by exchanging a photon 
(with a small cross section, the exchanged parti- 
cle might be also a W^, Z gauge boson). The 
kinematics of the reaction can be characterized 
by the following quantities 



• photon virtuality: Q 2 



-(k-k') 



ty2 



• Bjorken's scaling variable: x = Q 2 /(2P ■ q) 

• inelasticity: y = P ■ q/(P ■ k). 



2 Thereafter the term 'lepton' is used to refer either to 
electrons or positrons accelerated at HERA. 



previously explored kinematical range by about 
two orders in magnitude in x and Q 2 . 

3. Structure function and low-ir. parton 
densities 

The structure functions Fi and Fl are defined 
by the differential lepton-proton cross section 



' ep 



dxdQ 2 



(l + (l-y) 2 )F 2 (x,Q 2 ) 



Q 4 x 
-y 2 F L ( Xl Q 2 ) 



(1) 



where a denotes the fine structure constant. For 
most applications the longitudinal structure func- 
tion Fl can be neglected since y is sufficiently 
small. 
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In leading order perturbative QCD the struc- 
ture function F2 is related to the quark densities 
q, q in the proton by 



(qi(x,Q 2 ) + qi(x,Q 2 )), (2) 



where is fractional charge of the quark 
i. The Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) equations predict the Q 2 

dependence of the quark and gluon densities (for 
large Q 2 and not too small x) 



d 



dlnQ 2 



q( Xl Q 2 ) 
g{x,Q 2 ) 



2tt 



P 



q(x,Q 2 ) , ,,,, 



Here Py represents the splitting kernel describ- 
ing the probability of finding a daughter par- 
ton i in the parent parton j (e.g. J^ fe P(J — > 
i,k)) and we have used a s = g 2 /(4ir) with 
g s being the strong coupling constant. In 



the limit ln{l/x), ln(Q 2 /A 2 ) 



00 (double 



leading-logarithmic approximation) the DGLAP 
equations predict a steeply rising gluon density 
g(x,Q 2 ) 



xg(x, 



exp 
1 
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In- 



In 1 
-^ln- 



3 n f In $ 



(4) 



(5) 



where A denotes the QCD renormatization scale 
and nf is the number of quark flavours. Since 
the gluon and sea quark densities are closely re- 
lated (see Eq. (§)), a rise of F 2 is expected for 
decreasing x. This has been confirmed by HI 
and ZEUS measurements. Experimentally, the 
gluon density can be estimated, for example, from 
the scaling violation of F2 {Q 2 dependence) or 
from the diffractive J/^> production cross sec- 
tion. In Fig. the gluon density is shown for 
Q 2 = 20 GeV^ together with several theoreti- 
cal predictions. Currently, data on quark and 
gluon densities are available for x as low as 10 -4 
and Q 2 ~ 20 GeV 2 . These measurements lead 
to an improved knowledge of the low-x parton 
densities which has already been applied in many 
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NLO QCD fit 



CT = 20 GeV 

f I HI 1995+96 (0.118±0.005 

' ' (preliminary) m [ =1.3 - 1.8 Ge\ 

I I ZEUS 1994 (0.113±0.005 

' ' (preliminary) m.-1.3 - 1.5 Ge\ 

NMC (0.113) 



— MRSR1 (0.113) 
CTEQ4M (0.116) 
••• GRV94-HO (0.110) 




Figure 3. Gluon density at low x and Q 2 = 20 
GeV 2 . Shown are HI, ZEUS and NMC data to- 
gether with several parton density parametriza- 
tions. 



CR calculations, for example, calculations of neu- 
trino interaction length, charm production, and 
prompt muon yields ]T^ , [T^ | . 

Concerning ultra-high energy CR interactions, 
it is needed to extrapolate the measured parton 
densities down to lower x by several orders of 
magnitude. For example, partons with x values 
of 10~ 6 to 10~ 7 dominate minijet and charm pro- 
duction in proton-air interactions at £i a b ~ 10 20 
eV and are important in the case of neutrino-air 
interactions. The production of jets with several 
GeV transverse momentum in very forward di- 
rection involves even smaller x values of about 
HT 10 . 

From the HERA data it is clear that the gluon 
density rises rapidly at low x. It is convenient to 
parametrize this in terms of the power Ah with 
(see Eq. (I)) 



xg(x,Q 2 ) 



1 



(6) 



where data suggest Ajy ss 0.25 . . . 0.4. However, 
it is clear from simple geometrical arguments that 
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the rapid growth of the gluon density at low x will 
eventually be tamed by saturation effects [jl5],[l6| . 

Let's consider the scattering of a quark with 
the virtuality Q 2 off a proton. This virtual quark 
probes an effective transverse area of the size of 
~ a 2 {Q 2 )/Q 2 - As soon as the average trans- 
verse distance between the gluons in the proton 
becomes smaller than ~ f / y/Q 2 , several gluons 
participate in a single quark-proton interaction 
and the effective number of gluons " seen" by the 
quark is smaller than naively expected. At low x 
the effective number of gluons saturates with an 
upper limit of 



xg(x,Q 2 )~R 2 Q 2 



(7) 



where Rq is a measure of the radius of the trans- 
verse area available to gluons in the proton. This 
could be either the proton radius (i? 2 , ~ 5 GeV~ 2 ) 
itself or the size of gluon clouds around the va- 
lence quarks (hot spot scenario, R 2 ~ 2 GeV~ 2 ), 
see Fig. |[ It is important to note that satura- 



democratic distribution: 





hot spots: 





In 1/x ~ In s 



Figure 4. Two possible scenarios for the evolution 
of the parton density for decreasing x. 



simple geometrical considerations (see, for exam- 
ple, 1 19-^|; an overview on high-density parton 
evolution equations can be found in p^j). 

The gluon density is closely linked to the scal- 
ing violation of Fi 



xg(x,Q 2 ) 



4tt d.F 2 (f,Q 2 ) 
(40/27)a s dlnQ 2 



(8) 



This gives a handle to measure saturation effects. 
It follows from Eq. (Q) that in case of saturation 



dF 2 
d\nQ 2 



2^)2 



R Q 



(9) 



instead of dF2/d\nQ 2 



In Fig. 



ZEUS data are shown together with the predic- 
tion according to the DGLAP evolution equations 
(the Gliick-Reya-Vogt (GRV) f|| parton densi- 
ties have been used). Although it is possible to 
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tion effects might become important already for 
xg(x,Q 2 ) <C R 2 Q 2 as diffractive processes ob- 
served in DIS indicate (for example, see p7| , ^8[ ). 
Indeed, there are theoretical arguments support- 
ing the hypothesis that the gluon density might 
not reach the saturation limit (R) implied by such 



Figure 5. Fi scaling violation plotted as function 
of Bjorken-x. ZEUS data (circles) |24| are com- 
pared to the GRV prediction J|3| (triangles). 

fit the data within the DGLAP framework assum- 
ing a valence-like gluon density at low x p5|p6[ , 
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it is very likely that first hints of parton density 
saturation have been found [ p7| . 

The steep rise and the expected saturation of 
the gluon density at low x have several implica- 
tions for the interpretation of CR interactions. 
On one hand it is clear that the minijet cross sec- 
tion in hadron-air interactions increases rapidly 
with ^/s. On the other hand the saturation ef- 
fects make it very difficult to derive reliable cross 
sections for all processes where gluon densities at 
very low x are involved. This applies in partic- 
ular to the minijet and charm production cross 
sections in proton-proton and hence in hadron- 
air collisions and also, but to a lesser extent, to 
ultra- high energy neutrino cross sections. 

In most of the hadronic air shower simulation 
programs the production of minijets is one impor- 
tant ingredient. These Monte Carlo event genera- 
tors are based on the DGLAP evolution equations 
(|). Therefore it is interesting to consider the im- 
plication of the low- a; HERA data for the high- 
energy extrapolations done within these models. 
Without referring to some particular model, the 
HERA data imply that 

• the charged particle multiplicity rises with 
the energy faster than ln(s) but slower than 
s Ah , and 

• the mean transverse momentum of hadronic 
secondaries rises faster than ln(s). 

This can be understood by considering the 
transverse momentum distribution of hadronic 
secondaries produced in pp interactions as shown 
m Fig. | The region labeled " hard" is assumed to 
be calculable on the basis of the QCD-improved 
parton model. The "soft" region is dominated 
by saturation effects and cannot be calculated 
perturbatively at all. The transition region is 
effected by saturation effects but might still be 
calculable with perturbative means. For hard 
processes the differential jet cross section follows 
from 



da 



!(p5 



drjdp 2 ^ 



pi 



xig{xi,p\) x 2 g(x 2 ,p 2 j_) , (10) 



where x± t x 2 are the momentum fractions of the 
partons engaged in the hard scattering and p± is 




Figure 6. Transverse momentum distribution of 
hadrons produced in the central pseudorapidity 
region in CMS. 



the jet transverse momentum and r\ denotes the 
pseudorapidity. After integration over the phase 
space of the final state particles one finds, see for 



example 28 



> fco) 



(11) 



for a singular gluon distribution (^). Conse- 
quently, the hard part of the cross section rises 
much faster with the energy than the soft part 
(csoft ~ s As , A s = 0.07. ..0.15). In order to 
maintain a steady transition between soft and 
hard processes the value of fco bas to increase with 
the energy (which is also obvious from Eq. (0)). 
A simple estimate for the parameter ko can be 
made using pi 



fco ~ exp 



■fcVlns | , 



(12) 



which follows from Eqs. (^,0). The parameter c 
depends on the assumed low-as behaviour of the 
gluon density (^J). In Fig. [7] the results for c = 0.9 
(Ah = 0.3) and c = 1.12 {A H = 0.44) are shown. 

In a DGLAP-based model, realizing steeply ris- 
ing parton densities at low x but neglecting sat- 
uration effects, the average multiplicity of secon- 
daries becomes proportional to the minijet rate at 
very high energies. Consequently the multiplicity 
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log 10 (E lab /eV) log 10 (E lab /eV) 
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Figure 7. Estimate for the parameter ko charac- Figure 8. Average charged particle multiplicity 

terizing the transverse momentum below which in pp collisions. The data are taken from ]30|-|35[] 

saturation effects should be important. 



increases with s according to a power law 



rich ~ "jet 



CTjet 



„0.1...0.3 



(13) 



The inclusion of saturation effects reduces signif- 
icantly this rapid increase 



n c h 



k 2 



exp 



> ln(s) (14) 



but the multiplicity still grows considerably faster 
than ln(s) as implied by limiting fragmentation 
and Feynman scaling. Hence, HERA data ex- 
clude multiplicity extrapolations of the type (|l3| ) 
as found, for example, in QGSjet j29| but also 
disfavour a pure ln(s) extrapolation. In Fig. ||the 
two limiting extrapolations are shown, the upper 
curve corresponds to a DGLAP implied power law 
and the lower curve represents a ln(s) extrapola- 
tion of the collider data. 

In the very high energy limit the mean trans- 
verse momentum of secondaries is of the order of 
ko. Neglecting saturation effects (which means 
assuming fco =const.), the mean transverse mo- 
mentum of secondaries does not increase beyond 
a certain value even at ultra high energies. By 
contrast, the inclusion of saturation effects leads 
to an increase of the average p± which is faster 
than ln(s). The size of this effect can be seen in 
Fig. [|. Most of the currently used models show a 
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Figure 9. Average transverse momentum of sec- 
ondaries in pp collisions. Shown are model ex- 
trapolations of collider data 34 3^] with and 
without taking into account saturation effects. 
The calculation has been done with Phojet 



very moderate increase of the average p±_ which 
is comparable with the dashed curve J38[ . Only 
the Dpmjet Monte Carlo ]39| predicts a rapid in- 
crease of the average transverse momentum sim- 
ilar to the solid curve in Fig. ^[ In air shower 
experiments, the lateral spread of muons with en- 
ergies of the order of 1 GeV might allow to dis- 
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tinguish between the different model predictions. 

4. Hadronic final state in DIS 

Not only structure function data but also the 
knowledge gained at HERA on the hadronic final 
state of photon-proton interactions are of inter- 
est for CR physics. For example, hard interaction 
processes are governed by the same parton radi- 
ation processes in proton-air interactions and in 
DIS. As can be seen by comparing Fig. |l0| with 
f[l[ the parton emissions giving rise to the gluon 
which couples to the quark in the hard interac- 
tion are similar. The emission of partons (the 
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Figure 11. QCD-improved parton model view on 
jet production in proton-proton scattering. 
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Figure 10. QCD-improved parton model view on 
DIS. 



gluons shown in Figs. and |Tl|) is assumed to 
obey the DGLAP equations. Since basically al- 
most all perturbative calculations including the 
event generators used in CR physics are based on 
these evolution equations, an experimental confir- 
mation is very important]^] This is in particular 
true since there exist theoretical arguments that 
there might be a deviation from the DGLAP pre- 
dictions for low x processes po|-p^| . In the case of 
DIS the kinematics of the quark loop can be fixed 
by selecting events with certain scattering angles. 
By means of measuring the transverse momenta 

3 The description of the Q 2 -scaling violation of the struc- 
ture function F2 is an experimental confirmation of the 
DGLAP equations, however due to the unknown initial 
conditions not fully conclusive. 



of the gluon- induced jets fcr.i ■ ■ ■ fcr,3 the predic- 
tions of (H) can be tested. The DGLAP evolution 
equations imply 

xi > x 2 » x 3 . (15) 

In Fig. |l2| the results of the HI measurement 
for forward going jets are shown. Forward-jets 
roughly correspond to jets produced by the gluon 
with kx.i in Fig. |l0|. The inclusive cross section 
for these jets with 0.5 < k^/Q 2 < 2 is plotted 
as function of the Bjorken-x characterizing the 
DIS process. Clearly the DGLAP based calcula- 
tions cannot describe the data. Further measure- 
ments arc needed to confirm this deviation. There 
are some important consequences if the DGLAP 
based approximations break down at low x. 

• The incoming partons of hard scattering 
processes have sizable transverse momenta 
which cannot be neglected. 

• Strong kr ordering (|l5|) does not hold at 
low x. 

• One parton emission chain can give rise to 
more than two minijets having almost the 
same transverse momentum. 

The last conclusion can be mathematically re- 
phrased to 

1 f ^^"ict j9 / \ 
Chard ^ r / -J 2 d P_L,jet , (16) 

2 J d^Pi.jet 



(xlO" 2 ) 



Figure 12. Comparison of HI data on forward-jets 
with theoretical predictions. The dashed and dot- 
ted curves correspond to DGLAP based Monte 
Carlo calculations. The solid curve and the stars 
represent alternative approaches. For a detailed 
discussion see M. 



where Chard is the cross section for hard processes 
and daj C t/d 2 p±j C t is the inclusive jet cross sec- 
tion. Up to now all event generators applied in 
CR simulations assume that the hard cross sec- 
tion (which is needed for the calculation of the 
total cross section) is half the inclusive single- 
jet cross section given by perturbative QCD and 
collinear factorization (DGLAP approach). The 
size of the effect has been estimated in jh| and is 
shown in Fig. [l3| 

Concerning minimum bias physics, photon- 
proton collisions are very similar to pion-proton 
or proton-proton interactions. This is well under- 
stood in terms of the Vector Dominance Model 
(see, for example, |Q). The transition form real 
(Q 2 « 0) to highly virtual photons changes the 
event properties. The hadronic final state can 
be subdivided into photon fragmentation, cen- 
tral, and proton fragmentation regions. As has 
been confirmed experimentally, hadron produc- 
tion in the proton fragmentation region is almost 



Figure 13. Comparison of the hard cross section 
as calculated within the collinear factorization ap- 
proximation (dashed line) and an alternative ap- 
proach which is expected to be more suited for 
low- a; processes (solid line). The graph is taken 
from H. 



independent of the photon virtuality Q 2 . This 



supports the hypothesis of factorization |45| and 
allows us to compare data on leading baryon pro- 
duction in jp collisions to Monte Carlo predic- 
tions for pp interactions at the same energy. In 
Fig. [14] the energy fraction x lab = E haxyon /E heam 
carried by the leading proton in photon-proton 
and proton-proton collisions is shown . There 
is no Q 2 -dependence of the data within the exper- 
imental errors. The data show that there is no dip 
close to a:i a b ~ 0.9 as might be expected from the 
naive application of Feynman scaling violation ar- 
guments. The Sibyll HH] and QGSjet m 
models describe the data reasonably well whereas 
Dpmjet IH seems to underestimate the data 
close to iriab ~ 0.9. 

5. Summary and outlook 

HERA data and theoretical work related to 
HERA physics clearly improved our understand- 
ing of hadron production in such a way that the 
uncertainties in the model predictions for very 
high energy cosmic ray interactions can be re- 
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Figure 14. Energy fraction xi a b carried by the 
leading proton. The data are photoproduction 
and DIS measurements whereas the curves are 
model predictions for pp collisions. 



duced. 

Experimental evidence has been found for a 
steeply rising gluon density in the proton at low 
x. The inevitable consequence of this finding is 
that saturation effects are important. First ex- 
perimental signs for this have been reported. The 
implication of these results to high-energy extrap- 
olations done within models for CR interactions 
include 

• the particle multiplicity grows faster than 
ln(s) but slower than s Ah 

• the average transverse momentum of secon- 
daries rises faster than ln(s) 

Investigations of the hadronic final state in DIS 
at low x reveal deviations from the prediction 
obtained in the "standard" DGLAP framework 
(collinear factorization). Up to now, the physics 
of hard processes at low x is theoretically not well 
understood. However, it can be expected that 
the research activity on this subject triggered by 
HERA measurements will soon lead to a consid- 
erably improved understanding of typical low- a; 
processes such as minijet and charm production. 

It has been shown that HERA data can be used 
to reduce the uncertainties in Monte Carlo model 



predictions concerning central hadron or jet pro- 
duction as well as leading baryon distributions. 
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